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Formation of the Galactic globular clusters with He-rich 
stars in low-mass halos virialized at high redshift 
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ABSTRACT 

Recent observations have reported that the Galactic globular clusters (GCs) with 
unusually extended horizontal-branch (EHB) morphologies show a significantly lower 
velocity dispersion compared with that of the entire Galactic GC system. We consider 
that the observed distinctive kinematics of GCs with EHB has valuable information on 
the formation epochs of GCs and accordingly discuss this observational result based on 
cosmological N-body simulations with a model of GC formation. We assume that GCs 
in galaxies were initially formed in low-mass halos at high redshifts and we investigate 
final kinematics of GCs in their host halos at z — 0. We find that GCs formed in halos 
virialized at z > 10 show lower velocity dispersions on average than those formed at 
z > 6 for halos with GCs at z = 0. We thus suggest that the origin of the observed 
lower velocity dispersion for the Galactic GCs with EHBs is closely associated with 
earlier formation epochs (z > 10) of halos initially hosting the GCs in the course of the 
Galaxy formation. Considering that the origin of EHBs can be due to the presence of 
helium-enhanced second-generation stars in GCs, we discuss the longstanding second 
parameter problem of GCs in the context of different degrees of chemical pollution in 
GC-forming gas clouds within low-mass halos virialized at different redshifts. 
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1 INTRODUCTION 

A growing number of observational studies have recently 
reported that some Galactic GCs (e.g., u> Cen and NGC 
2808) show a possible He abundance spread in their stellar 
populations (e.g., Bedin et al. 2004; D'Antona et al. 2005; 
Lee et al. 2007, L07; Piotto et al. 2007). L07 furthermore 
have reported that about 25% of Galactic GCs show very 
extended horizontal-branch (HB) morphologies that can be 
due to He-rich stars in these GCs. One intriguing observa- 
tion of the Galactic GCs with extended HB morphologies 
(EHB) is that they show a significantly lower velocity dis- 
persion (93 ± 13 km s^ 1 ) in comparison with the rest of the 
halo GC system (137±14 km s" 1 , L07). Although this obser- 
vation can shed new light on the Galaxy formation through 
hierarchical merging of low-mass building blocks (L07), the 
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origin of the unique kinematics of the Galactic GCs with 
possible He-rich stars remains unclear. 

Although previous one-zone chemical evolution models 
have discussed the observed possible star-to-star variation in 
He abundance in Galactic GCs (e.g., D'Antona et al. 2002; 
Bekki et al. 2007a) , they have not discussed the origin of the 
observed kinematics of the GCs with EHB. Recent numeri- 
cal simulations of GC formation during hierarchical galaxy 
formation based on a ACDM model have greatly improved 
the predictability of dynamical and chemical properties of 
GC systems in galaxies (Kravtsov & Gnedin 2005; Yahagi 
& Bekki 2005; Bekki et al. 2007b). However they have not 
yet discussed the origin of the observed lower velocity dis- 
persion of Galactic GCs with EHB (hereafter these GCs are 
referred to as EHB GCs for convenience). 

The purpose of this paper is to demonstrate, for the first 
time, that the observed kinematics of EHB GCs is consistent 
with a scenario in which they were formed in low-mass halos 
virialized before z > 10. We consider that (i) the present-day 
metal-poor GCs of galaxies were formed initially in low-mass 



2 K. Bekki, H. Yahagi, M. Nagashima, and D. A. Forbes 



halos virialized at z > 6 before the completion of reioniza- 
tion (Fan et al. 2003) and (ii) GCs were then tidally stripped 
from their host halos during hierarchical merging to become 
halo GCs in the final galaxy. The above (i) is consistent 
with recent hydrodynamical simulations of GC formation in 
low-mass halos by Bromm & Clarke (2002). Based on high- 
resolution cosmological N-body simulations with a model of 
GC formation, we investigate the kinematics of GC subpop- 
ulations originating from low-mass halos virialized at differ- 
ent redshifts (z). We show that the kinematical properties 
of GC subpopulations initially in low-mass halos virialized 
well before z > 10 show significantly lower velocity disper- 
sions in comparison with GCs formed at z > 6. We discuss 
this result in terms of the origin of EHB GCs in the Galaxy. 



2 THE MODEL 

Since the present model of GC formation in hierarchical 
galaxy formation is essentially the same as those adopted in 
our previous works (Yahagi & Bekki 2005; Bekki et al. 2006; 
Bekki et al. 2007b), we briefly explain the model here. We 
simulate the large scale structure of GCs in a ACDM Uni- 
verse with n = 0.3, A = 0.7, H = 70 km s _1 Mpc -1 , and 
(jg = 0.9 by using the Adaptive Mesh Refinement A— body 
code developed by Yahagi (2005) and Yahagi et al. (2004), 
which is a vectorized and parallelized version of the code de- 
scribed in Yahagi & Yoshii (2001). We use 512 3 collisionless 
dark matter (DM) particles in a simulation with the box 
size of 70fr~ 1 Mpc and the total mass of 4.08 x 1O 16 M . We 
start simulations at z — 41 and follow it till z — in order 
to investigate physical properties of old GCs within virial- 
ized dark matter halos at z = 0. We used the COSMICS 
(Cosmological Initial Conditions and Microwave Anisotropy 
Codes) , which is a package of fortran programs for generat- 
ing Gaussian random initial conditions for nonlinear struc- 
ture formation simulations (Bertschinger 1995, 2001). 

Our method of investigating GC properties is described 
as follows. Firstly we select virialized dark matter subhalos 
at each output by using the friends-of-friends (FoF) algo- 
rithm (Davis et al. 1985) with a fixed linking length of 0.2 
times the mean DM particle separation. The minimum parti- 
cle number AT m i n for halos is set to be 10. For each individual 
virialized subhalo, the central particle is labeled as a "GC" 
particle. This procedure for defining GC particles is based 
on the assumption that energy dissipation via radiative cool- 
ing allows baryon to fall into the deepest potential well of 
dark-matter halos and finally to be converted into GCs. We 
assume that GC formation is truncated for halos virialized 
after z = z t and that z% = 6 corresponds to the epoch of 
the completion of reionization (e.g., Fan et al. 2003). Sec- 
ondly, we follow the dynamical evolution of GC particles till 
z = and thereby derive locations (x,y,z) and velocities 
(v x ,Vy,v z ) of GCs at z — 0. We then identify virialized ha- 
los at z = with the FoF algorithm and investigate whether 
each GC is within a halo. 

We investigate velocity dispersions (er Zv<z ) of GC sub- 
populations originating from halos virialized at redshifts 
z v < z in halos at z = 0. The most important key pa- 
rameter in the present study is the virialization redshift z v , 
and accordingly we investigate the kinematics of GC sub- 
populations with 6 < 2 V < 15 for z% = 6 in the GC system 



(GCS) in each of the simulated halos at z — 0. We partic- 
ularly investigate the ratio of a Zv<z of GC subpopulations 
to the entire GCS in each halo at z = 0, denoted as (7 a ii- 
These velocity dispersion ratios are referred simply to as R B 
from now on. To avoid a large error bar of Rs (=a Zv<z /a a _n 
for each individual halo at z = 0) resulting from the small 
GC number, we pick up haloes in which both the numbers 
of GCs originating from halos with z v < z and z% < z < z v 
exceed a threshold number of n th . When the total number 
of GCs in the i-th halo is n gc ,i, this criterion corresponds to 
n gCj i > 2n t h- We mainly show the results for the models with 
nth = 3 in which n gCji in each halo is larger than 6 (= 2n t h)- 
In order to discuss the dependences of the present results on 
nth, we also show the results of the models with nth = 24. 
Although an error bar in Rs for each individual halo be- 
comes significantly small in the models with nth = 24, the 
total number of halos with nth > 24 becomes very small. 

L07 reported that the Galactic GCs with EHB has a 
velocity dispersion of 93 ± 13 km s - (= ctehbgCs) whereas 
the entire GCS of the Galaxy has a dispersion of 124± 10 km 
s _1 (= a a ii). We compare the observed ratio (R ) of oehbgcs 
to a a ii (shown in Table 1 of L07) for the Galactic GCS with 
our simulated one in order to determine the best value of 
z v for which the observed R Q can be well reproduced. Since 
the observed <T a ii is derived for all GCs including relatively 
young halo GCs (L07), z% = 6 is regarded as a reasonable 
truncation epoch in deriving a a ii for the present study. 

The Galactic GCs around the Galactic bulge ("bulge 
GCs" or "disk GCs") could have been formed as a re- 
sult of dissipative (and star-forming) major merging pos- 
sibly responsible for the bulge formation at a lower redshift. 
The present model does not include the formation of these 
younger, metal-rich GCs, and accordingly a&n is estimated 
only for GCs originating from halos virialized before z~6. 
Thus it should be stressed that the derived cr a ii in the present 
model does not literally mean the velocity dispersion of all 
GCs in a galaxy. However we think that it is still reasonable 
to compare the derived (7 a ii in the simulations with the ob- 
served aaii (based on all GCs with known HB morphologies), 
because the observed small number fractions of bulge/disk 
GCs in disk galaxies (e.g., Zinn 1985 for the Galaxy) mean 
that there could be only a small difference in velocity dis- 
persions between a GCS including bulge/disk GCs and the 
GCS not including them for a galaxy. 



3 RESULTS 

Fig. 1 shows how Rs for GC subpopulations with z v = 7, 10, 
and 13 depend on Mh in the model with nth = 3 correspond- 
ing to n gCj i > 6. Fig.l demonstrates that GC subpopulations 
with higher z v are more likely to show a lower R s for a given 
Mh, though the dispersion in R B for a given Mh is larger for 
low-mass halos. This result suggests that R s of a GC sub- 
population in a galaxy can provide some information about 
the redshifts of virialization for low-mass halos that initially 
hosted the subpopulation. Fig.l also shows that more mas- 
sive halos are more likely to show a higher Rs value for a 
given z v , which means that differences in velocity dispersions 
between GC subpopulations originating from high-z halos 
with different redshifts of virialization are less remarkable 
in more massive galaxies. 
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Figure 1. Velocity dispersion ratios (R s = cr Zv < z /cr a ii) as a 
function of their host halo masses (M^) for GC subpopulations 
with the virialization redshifts z v = 13 (top), 10 (middle), and 
7 (bottom). Here the ratios are derived for the ^-components 
of velocity dispersion (<r x ) in halos at z = 0. Each small dot 
represents a halo with a GCS with n gCi ; > 6. Blue, filled circles 
represent the mean values of the ratios in five halo mass bins and 
the error bars are due to the number of halos (N^) in each halo 

mass bin (i.e., oc 1/v/ 2(N\ l — 1)). For the lowest mass bin with 
no GCs for z v = 13, o- Zv < z /cr a ii=0 is shown. For comparison, 
the observed ratio (Ro) of o"EHBGCs/o" a ii (^07) is shown by a 
red, open circle for the mass model of the Galaxy (Wilkinson & 
Evans 1999). Note that GC subpopulations with z v = 10 and 
Mb ~ 10 12 Mq show low o" Zv < z /cr a ii which is best agreement to 
the observed one in L07. 



Furthermore, Fig. 1 shows that GC subpopulations with 
z v = 10 in galaxy-scale halos with Mh w 1O 12 M0 (cor- 
responding roughly to the total mass of the Galaxy, e.g., 
Wilkinson & Evans 1999) have _R S very similar to the ob- 
served velocity dispersion ratio (R ): The Rs of GC sub- 
populations with z v — 7 are too high to be consistent with 
R , whereas those of GC subpopulations with z v — 13 are 
marginally consistent with R a . These results suggest that 
EHB GCs in the Galaxy were initially within low-mass ha- 
los virialized at z > 10 and later tidally stripped during the 
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Figure 2. The same as Figure 1 but for <r z (top), cr y (middle), 
and (T x (bottom) of GC subpopulations with z v = 10. 



1.5 



b 0.5 



1 

"«... >SM 


1 1 1 1 1 1 

• • 


1 c 

1 1 1 1 ill 


1 1 1 1 1 1 1 1 1 



11 12 13 

log 10 M h (M ) 



14 



Figure 3. The same as Figure 1 but for GC subpopulations with 
z v = 10 in the model with n gc ; > 24. 
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hierarchical merging of halos to finally become the halo GCs 
in the Galaxy. The origin of the lower R s in GC subpopu- 
lations with z v = 10 may well be closely associated with 
characteristic orbital properties (e.g, more eccentric orbits 
and smaller pericenter distances with respect to their halos' 
centers). 

Fig. 2 shows that lower R B values for GC subpopula- 
tions with z v — 10 can be clearly seen in the three velocity 
components (i.e., a x , a y , and <r z ) in galaxy-scale halos with 
Mh ~ 1O 12 M0, which confirms that lower R B values can be 
one of the kinematical properties characteristic of GC sub- 
populations with z v = 10 in galaxy-scale halos. Fig.2 also 
shows that the trend of GC subpopulations in more mas- 
sive halos (10 Mq < Mh) to have higher R s values can be 
clearly seen for the three velocity components. This result 
implies that differences of line-of-sight velocity dispersions 
in GC subpopulations (e.g., intra-group and intra-cluster 
GCs) originating from halos virialized at different redshifts 
are observationally difficult to detect for group-scale and 
cluster-scale halos. 

Fig. 3 shows R s of halos in the model with n th = 12 cor- 
responding n gc ,i > 24 in which only halos with significantly 
larger numbers of GCs with z v = 10 are selected. Lower R s 
values in halos with Mh w 1O 12 M0 can be clearly seen even 
in this strongly biased sample of halos at z — 0, though the 
scatter in R B become small thanks to a larger number of 
GCs used in deriving _R S . This result demonstrates that the 
derived lower R s values in galaxy-scale halos are not due to 
the parameter nth introduced for minimizing scatters in R B 
for individual halos in the present study. 

4 DISCUSSION AND CONCLUSIONS 

4.1 Origin of GCs with He-rich stars 

Recent observations have revealed that the (line-of-sight) ve- 
locity dispersions of "ultra-compact dwarf" (UCD) galaxies 
are significantly smaller than those of other galaxy popula- 
tions in the Fornax and the Virgo clusters of galaxies (Mieske 
ct al. 2004; Jones et al. 2006; Gregg et al. 2007). Bekki (2007) 
first demonstrated that the observed lower velocity disper- 
sion of the UCD population in the Fornax cluster is consis- 
tent with the UCDs having significantly smaller pericenter 
distances (with respect to the center of the Fornax) than 
other galaxy populations in the cluster. This result suggests 
that the simulated lower velocity dispersions of GCs with 
z v > 10 in halos at z = are due to significantly smaller 
pericenter distances of the GCs (with respect to the halos' 
centers) in the present study, ft furthermore implies that 
EHB GCs in the Galaxy can have smaller pericenter dis- 
tances with respect to the Galactic center. 

The present study has first shown that EHB GCs origi- 
nate preferentially from low-mass halos that were virialized 
at higher redshifts (z > 10). Given the fact that the origin 
of EHB can be closely associated with the presence of He- 
rich second generation stars in GCs (e.g., D'Antona et al. 
2002; L07; but see Choi & Yi 2007 for a different formation 
mechanism of He-rich stars), the above result implies that 
the formation of He-rich stars is much more likely to happen 
in GC-forming gas clouds of halos with possibly higher gas 
densities (due to higher z v ). Gas ejected from the first gener- 
ation of stars of a GC can be efficiently trapped in the central 



regions of halos during the GC formation owing to the deep 
gravitational potential of the halos so that the gas can be 
used for the formation of the second generation with a higher 
helium abundance for the GC (Bekki 2006): there are two 
generations of stars with different helium abundances in a 
GC. Formation of He-rich second-generation stars therefore 
can proceed more efficiently in the central regions of halos 
with higher z v and thus deeper gravitational potentials (for 
a given halo mass). Therefore, formation of He- rich stars is 
highly likely to occur in GC-forming gas clouds located in 
the central regions of halos with higher z v . 

Although both observational and theoretical studies 
have discussed the origin of the second parameter(s) that 
can control the morphologies of HBs in GCs, it is not yet 
clear what the second parameter(s) is (are) in these stud- 
ies (e.g., Catelan et al. 2001 and references therein). Re- 
cently, Recio-Blanco et al. (2006, R06) have discovered that 
HB morphologies depend on the total masses of GCs in the 
sense that more massive GCs tend to have more extended 
HBs: the masses of GCs can be one of the second parame- 
ters. They have accordingly suggested that higher fractions 
of He-rich stars in more massive GCs can cause EHB in the 
GCs and that the origin of the higher fractions can result 
from more efficient "self-pollution" by stars formed earlier 
in more massive gas clouds. 

If the number fraction of He-rich stars in GCs is one 
of the second parameters, as suggested by R06, then the 
present numerical study suggests that the differences in HB 
morphologies in GCs are due partly to the differences in 
formation epochs (i.e., the redshifts of virialization) of low- 
mass halos initially hosting GCs. The present study further- 
more suggests that if more massive GCs with He-rich stars 
are more likely to be formed in halos virialized at higher red- 
shifts, then the observational results by R06 (i.e., the pres- 
ence of EHB in massive GCs) and by L07 (i.e., kinematic 
decoupling of EHB GCs) can be self-consistently explained. 
Owing to the lack of extensive theoretical studies on GC 
formation processes in halos with different virialization red- 
shifts, it is however unclear whether more massive GCs can 
be formed in halos virialized at higher redshifts. 

It remains unclear whether AGB stars (e.g., D'Antona 
ct al. 2002; Karakas et al. 2006; Bekki et al. 2007) or mas- 
sive stars (e.g., Charbonnel & Prantzos 2006; Bekki & Chiba 
2007) are responsible for the origin of He-rich stars in GCs. 
Miocchi (2007) has recently discovered an intriguing corre- 
lation between the possible presence of intermediate-mass 
black holes (IMBHs) and the presence of EHB in GCs. 
Given the fact that recent numerical simulations (e.g., Porte- 
gies Zwart et al. 2004) have demonstrated the formation 
of IMBHs through runaway collisions of massive stars, the 
above result by Miocchi (2007) implies that massive stars 
rather than AGB stars could have played a major role in 
the formation of He-rich stars within forming GCs. 

4.2 Chemo-kinematics of GC subpopulations and 
galaxy formation 

Although previous theoretical and observational studies 
have focused mostly on differences in kinematical properties 
between metal-poor and metal-rich GCs in galaxies (e.g., 
Bekki et al. 2005; Brodie & Strader 2006), they have not 
so far investigated the kinematical differences between GCs 
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with different masses, abundances of He and light elements 
(e.g., C, N, and O), and cv-abundance ratios (e.g. [Mg/Fe]). 
The discovery of the chemokinematical correlation in the 
Galactic GCs by L07, combined with the present numerical 
results (i.e., lower velocity dispersions of GC subpopulations 
originating from halos virialized at higher redshifts), pro- 
vide the following three implications. Firstly, there could be 
kinematical differences in GCs with He-rich stars compared 
to those without in many galaxies. Recently, Kaviraj et al. 
(2007) have investigated the ultraviolet and optical proper- 
ties of 38 GCs in M87 and suggested the possible presence 
of He-rich stars in the GCs. It is thus an interesting obser- 
vational question whether the 38 GCs in M87 show a lower 
velocity dispersion in comparison with the entire GC popu- 
lation in M87. 

Secondly, velocity dispersions of GC subpopulations 
with very high [a/Fe] (e.g., [Mg/Fe]) can be lower. This is 
mainly because star formation time scales, which determine 
[a/Fe] of stellar populations (e.g., Matteucci & Frangois 
1992; see also Nagashima et al. 2005), is shorter in halos 
that were formed at higher redshifts and thus have higher 
densities in their baryonic components. Thirdly, there can 
be significant kinematical differences between GC subpopu- 
lations with lower masses and those with higher masses in 
the sense that the latter subpopulations should show lower 
velocity dispersions. This is because the present simulations 
suggest that more massive GCs can be formed at halos viri- 
alized at higher redshifts. 

The observed chemokinematical correlations in the 
Galactic stellar halo and the physical properties of GCs 
have long been discussed in the context of the time scale of 
the Galaxy formation (e.g., Eggen, Lynden-Bell, & Sandage 
1962) and the accretion history of dwarfs with GCs (Searle 
& Zinn 1978) . The present study suggests that the observed 
kinematical correlations of GCs with, and without, EHBs 
in the Galaxy can be considered in terms of different for- 
mation epochs of low-mass halos initially hosting GCs. We 
thus suggest that HB morphologies of GCs in a galaxy con- 
tain valuable information both on the formation epochs of 
building blocks of the galaxy and on the hierarchical merg- 
ing histories of that galaxy. We also suggest that one of the 
second parameters governing the HB morphology is the for- 
mation (i.e., virialization) epochs of low-mass dark matter 
halos that initially hosted GCs at high redshifts. In other 
words, the environments of GCs are a key 'second parame- 
ter' as they determine the extent to which GC-forming gas 
clouds are chemically polluted by earlier generations of stars. 
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